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^ \ Abstract: We compute three-body nuclear forces at short distances by using the 

• \ nuclear matrix model of holographic QCD proposed in our previous paper with P. Yi. 

O \ We find that the three-body forces at short distances are repulsive for (a) aligned 



three neutrons with averaged spins, and (b) aligned proton-proton-neutron / proton- 
neutron-neutron. These indicate that in dense states of neutrons such as cores of 
neutron stars, or in Helium-3 / tritium nucleus, the repulsive forces are larger than 



^ ' the ones estimated from two-body forces only. 
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1. Introduction. 

One of the fundamental ingredients of nuclear physics is the nuclear force with which 
point-like nucleons interact with each other. A variety of aspects of nuclear forces 
results in the protean metamorphosis of nuclei, the bound states of nucleons. It 
is known that in nuclear forces there are forces that can not be explained by two- 
body forces only, one of which is the three-body force. The three-body forces play 
important role, for example, in reproducing excitation spectra of light nuclei, or 
explaining equations of states for high-density baryon matters such as supernovae 
and neutron stars. However, in spite of the long history of nuclear physics, the bulk 
properties of three-body nuclear forces are yet to be revealed. 

The main obstacle for revealing the various aspects of nuclear forces is obvious: 
QCD is strongly coupled and thus difficult to solve. In this paper, by using a nuclear 
matrix model of holographic QCD which we have derived together with P. Yi in [|^ , 
we explicitly compute a three-body nuclear force in a large Nc holographic QCD. 
The two-body nuclear force was already computed in [Q. 

For the derivation of our matrix model [0 we use the gauge/string duality (the 
AdS/CFT correspondence) 0, ^, Q applied to a D4-D8 system [^, ^ of a large N^ 
QCD at a large 'tHooft coupling A. Precisely speaking, our matrix model is a low- 
energy effective field theory on baryon vertex D4-branes in the D4-D8 holographic 



model [^, ^ of large Nc QCD. The matrix model describes /c-body baryon systems 
with arbitrary k, where the size of the matrix is given by this k, based on the fact that 
baryons are wrapped D-branes on sphere in the gravity side of the gauge/string 
duality. In the previous work |I|, in addition to the derivation of the matrix model, 
the cases with k = 1 (baryon spectrum) and k = 2 (two-body nuclear force) were 
studied. For fc = 2, it was found that a universal repulsive core exists for any baryon 
states with two flavors. Since our matrix model is not a phenomenological model 
for multi-baryon systems, but based on a firm ground of the gauge/string duality in 
string theory, it is natural to extend the analysis of our matrix model to derive the 
three-body nuclear forces. In this paper, we continue the analysis to the k = 3 case, 
i.e. we study the short-range three-body nuclear force, using the matrix model. 

Although generic configurations of three baryons can be treated in the matrix 
model, as the computations are involved and thus not so illuminating, in this paper we 
shall concentrate on two particular examples: (a) three neutrons with spins averaged 
and (b) proton-proton-neutron (and proton- neutron-neutron) , both aligned on a line 
with equal spacings. System with spin averaged is rather typical for dense states of 
multi-baryons such as cores of neutron stars. The latter is related to Helium-3 and 
tritium nucleus. For both cases, the resultant three-body potential is positive, i.e. 
repulsive. It scales as N^/X'^r^ (where r is the inter- nucleon distance), in contrast to 
the two-body repulsive core ~ Nc/Xr"^. As the region of validity is at short range, 
1/VXMkk < r < 1/Mkk (where Mkk ~ C(l GeV)), the three-body potential is 
suppressed compared to the two-body potential by ~ 1/A (rM^i^)^ <^ 1. However at 
very short distances, i.e. at high dense states of nucleons, three-body forces are not 
small. 

The organization of this paper is as follows. We first review the matrix model and 
the two-body calculation shown in |1|. Then in section 3, we calculate the three-body 
forces. First, as an exercise, we treat the case with spin/isospins aligned classically, 
and find that the three-body force vanishes for this case, which is consistent with 
the soliton approach P]. After that, we proceed to generic three-body forces with 
quantum spin/isospins. The final section is devoted for discussions. 



2. Review: a matrix model and two-body interactions. 

The procedures of the computation of the three-body nuclear forces is quite analogous 
to the two-body case performed in [|l|. Here we provide a summary of the matrix 
model action and the computation of the two-body nuclear forces of ^ . 

2.1 A matrix model action. 

In |T[, we proposed with P. Yi a U{k) matrix model which describes generic /c-body 
interaction of nucleons. Note that the rank of gauge group U{k) is not at all related 



to the number of colors N^. but just the number of nucleons k. The matrix model 
action is quite simple. 
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The peculiar property of this matrix model is the simplicity: changing the number 
of the nucleons k is available just by choosing U{k) for the gauge group of the matrix 
model. In [|T| it was demonstrated how to compute the baryon spectrum {k = 1) and 
two-body nuclear forces {k = 2) at short distances. 

To be concise, here we briefly describe the matter content of the matrix model 
( p.l|) . The model has a unique scale Mkk, and A = A^c5'qcd i^ ^^^ 'tHooft coupling 
constant of QCD, with the number of colors Nc- The fleld content is summarized in 
the following table. 
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The dynamical fields are only X*^ and Wai, while Aq and Dg are auxiliary fields. In 
writing these fields, the indices for the gauge group U{k) are implicit. In this paper 
we consider only the two-fiavor case Nf = 2 for simplicity. The symmetry of this 
matrix quantum mechanics is U{k)\ocai x SU{Nf) x S0{3) where the last factor S0{3) 
is the spatial rotation, which, together with a holographic dimension, forms a broken 
5*0 (4) ^ SU{2) X SU{2) shown in the table. The breaking is due to the mass terms 
for X'^ and w^i- In the action, the trace is over these U{k) indices, and the definition 
of the covariant derivatives is DqX^ = OqX^^ — i[AQ, X*^], Z^o""^ = dow — iwAo and 
Z^o""^ = dow + iAow. The spinor indices of X are defined as Xaa = X^^{aM)aa and 
j^aa ^ j^M^^^^^aa -y^^i^gre (Tm = (iT, 1) and (Jm = (— ir, 1), with Pauh matrices r. 
All of these definitions follow the notation of . For the derivation of this matrix 
model via gauge/string duality, see [Q. 

2.2 Two-body nuclear forces. 

We review briefiy [Q for explaining how to obtain the two-body nuclear forces. First, 
we describe a single baryon wave function, and second, obtain the two-body Hamil- 
tonian by integrating out Aq with a simple gauge choice. 

In all the cases, we need to solve the "ADHM constraint" iTOl which minimizes 



the potential induced by integrating out Dg. This is because the potential has a 



coefficient A^ which is very large in the gauge/string duahty. 
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Here A, B = l,...,k. 

For a single baryon k = 1, this equation is simply solved by Wai = UaiP where U 
is an SU{2) matrix and p is a constant. After integrating out the auxiliary field Aq, 
the matrix model action becomes a standard quantum mechanics whose Lagrangian 



is almost the same as that of the soliton approach |113| , [14 
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This quantum mechanics is solved by following |TB, Ti]. At the leading order in the 
large Nc limit, the wave functions for X'^ and p are classical, which simply means 
that we can substitute the classical values 

X^ = , p = 2-1/^3^/^ v^A-^/^Mk^ . (2.4) 

The wave functions for the spin/isospin U is nontrivial. They are shared with those 



of the Skyrme model |TT], |T2[, as described in ||T^, and given by 
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The SU{2) matrix U is represented by a unit 4-vector a as f/ = zajX* + a4l2x2j with 
(ai)2 + {a2f + {asf + {a,f = 1. 

Next, let us review the case /c = 2 of lU. The generic solution to the ADHM 
constraint (g^) with k = 2, Nf = 2 is the well-known ADHM data of two SU{2) YM 
instantons. 
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Here the off-diagonal part Y is defined as 
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The vector r (M = 1,2,3,4) is the distance between the two baryons, and \r\ = 
(r^^y. The SU{2) matrices U^^^ and f/*^^^ denote the moduli parameters of each 
baryon, while pi and p2 denote the moduli parameter associated with the size of 

— * 

instantons of each baryon. With this choice, the potential associated with D [i.e. 
the ADHM potential) vanishes. 



We integrate out the auxiliary field Aq to obtain the two-body Hamiltonian. 
With the U(2) decomposition Aq = A^12x2+AIt^+AIt'^ + AIt^, it is straightforward 
to evaluate the terms including the gauge field Aq in the action, 
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5ki„.+cs = ^^^^ [dt tr [{DoX^'f + DowtDoWo.,] + N, j dt tiAo 

2{Alf\r\' + 8{AinY'r + 2(p? + pl) {{<)' + K)^ + {Alf) 
+Ap,p,A',Al tr [f/«tf/(2)] + A{pl - pl)A',Al + -^A^l . (2.8) 

Solving the equations of motion for all the components of Aq and substituting the 
solution back to this action, we obtain the potential J dt V = — SkuL+cs; 

V = 2Vi_body + V2_body , ^1-body = TTTT n ' (^-Q) 

4AiWKK P 

l/2-body - ^^^^ |^|2 ^ 2p2 _ 2(^o)2p2 ' ^^'^^J 

Here Mq = (l/2)tr [[/(^)''"[/(^)] , and we put pi = P2(= p) which is justified as we keep 
only the leading term in the large N^ expansion. The value of p is ( |2.4| ).^ 

In addition to the terms in S'kin.+cs, there is the mass term for X'^ in the La- 
grangian, 

^^^ ■ iM^KtriXr = ^M|k {{rr/2 + 2(y^)^) . (2.12) 

The off-diagonal component Y gives an additional two-body potential. 
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■zr^(f/«)V(^)])', (2.13) 



where i = 1,2,3. So the 2-body potential is a sum of (|2.10|) and ( p.l3|) . The four- 
dimensional distance |r| is equal to the inter-baryon distance |r*p in three dimensions, 
since the classical value of the X^ for the single instantons is zero at the large Nc 
leading order, as in (p^). 

Using the nucleon wave function (|2.5| ), it is straightforward to evaluate the vac- 
uum expectation value of this potential. The final form of the two-body nuclear 



^For our later purpose, we write the expression for the case of classically aligned spins and 
isospins. This corresponds to t/'-^-' — U^^\ which is nothing but an ADHM data for 'tHooft 
instantons. The two-body potential is found as 



ci ^InNr, 1 



^'^-''""y ~ XM^W ^^'^^^ 



potential is (^)/i,ji,/2,J2 — ^c{^ + 5'i2Vt(^) with the standard definition S12 = 
12 J^fVgr-' — 4J{J2 (with f* = r^/\r\,i = 1,2,3), where the central and tensor forces 



are^ 



Mr-) = .('f H-8/;/W^^)^^. MO = 2./;/.^^. (2.14) 



3. Three-body baryon interaction. 

The three-body interaction potential can be computed by using the matrix model 
with k = 3 for k x k matrices. The procedures to compute the nuclear potential are 
parallel to the case of the two baryons in the previous section, and here is a summary 
of the procedures: 

(1) Choose your k, and solve the ADHM constraint ( p.2| ) (which minimizes the 
potential obtained by integrating out the auxiliary field D). 

(2) Substitute the solution back to the matrix model Lagrangian. 

(3) Integrate out the auxiliary field Aq. 

(4) Evaluate the Hamiltonian with your favorite baryon state. The baryon state 
is just a tensor product of single-baryon states (which is given by the k = 1 
analysis) . 

In this section, at first as an exercise, we consider a case where all the three baryons 
share the same classical spin/isospins. Then for next, after giving an explicit set-up 
for generic quantum spin/isospins, we demonstrate exact computations for baryons 
aligned on a straight line with equal spacings. The reason for choosing this linear 
position is just to simplify and illuminate the computations. Finally we evaluate 
the three-body Hamiltonian with specific three-baryon quantum states: (a) three 
neutrons with spins averaged, and (b) proton-proton-neutron and proton-neutron- 
neutron. We find that the three-body nuclear potential is positive i.e. repulsive. 

3.1 An exercise : classical treatment of spin/isospin. 

Let us evaluate the three-body Hamiltonian, first, for a simple situation where all 
the baryons share the same classical spin/isospin, as an exercise. 

Procedure (1): solving the ADHM constraint 

First, let us consider the configuration space of minimizing the D term. This is 
equivalent to the so-called ADHM constraints ( p.2|) for any A,B = 1,2,3. A simple 



^The result is quite close to that of the soliton approach [ha. 



solution to this constraint equation is the ADHM data for 'tHooft instantons, which 
we treat in this subsection, while in the later subsection we consider generic solution 
to this constraint. The ADHM data for the 'tHooft instantons consists of diagonal 
matrices X and special w's sharing the same orientation, 

wi = UMP^ (A = 1,2, 3), (3.1) 

where 2x2 unitary matrix U is independent of the index A. As the degrees of 
freedom w correspond to the spin and the isospin, this means that all the three 
baryons share the same "classical" spins and isospins. Here, fixing the orientation 
U for baryons cannot be achieved with wave functions with finite width, that is the 
reason we call this "classical." As the off-diagonal elements of X*^ vanish, all the 
commutators [X, X] are zero, which trivially satisfies (p.2|) . 

In p|, 'tHooft instantons are used in the soliton approach to evaluate the three- 
body nuclear forces. The result turns out to vanish. In this subsection, we will find 
that our matrix model also gives the same answer, the vanishing three-body force 
for the ADHM data of the 'tHooft instantons. 

We are going to choose implicitly the gauge d^w^^ = so that there is no time- 
dependence in w. See |T]] for details of the gauge choices. The matrices X^^ whose 
diagonal elements with M = 1,2,3 specify the spatial location of the baryons are 
diagonal, 

^*^=Ey^''- (3-2) 

a=3,8 

Procedure (2): substitute the ADHM data to the action 

The inter-baryon potential consists of two terms, the potential coming from the 
integration of Aq, and the potential from the mass term for X^. The latter vanishes 
for the 'tHooft instantons, as there is no off-diagonal extra component in ( |3.2| ). So 
in this subsection we concentrate on the former. 

Given the ADHM data, we can integrate out Aq, in analogy to the two-body 
case. The auxiliary field Aq is expanded by the Gell-Mann matrices A'^, 

Ao = A°l3x3 + E^oY- (3.3) 

a=l 

As in the two-body case, in the Lagrangian the terms containing Aq are (DqX)^ and 
DqwDw. The CS term contains only the overall U{1) component, Aq. 
First, the kinetic terms of X gives 

triDoX^r = \ {{Alr^r + (^^-3^^^) + ^ {«)' + (^^)^) (rf + ^^4^ 

+1 {(Air + (Air) {rr - V3r^r ■ (3.4) 



Next we consider the kinetic term for w, 

tiDowfDow^. = {A',Y (y1 I^^^I' ) + 2^S^o/ + i^'A^.Al (3.5) 

where a(= 1, ■ ■ ■ ,8) is the adjoint index of the SU{3), and A{= 1, 2, 3) is the index 
for the baryons. The coefficients j" and t"''' are defined as 

?-^f.^«.)% i'^'^^i^^^i^Zr- (3.6) 

Using the definition of symmetric structure constants d"''"^ for SU{3), 

the term with t"'' is replaced by 

1 Jabc 

i'^'A^X = -w^.iwi.yS'^'A^.A'^, + —wiX:,ciy^grA^o< ■ (3.8) 

Now due to the ADHM data (p.l|), all w^^ are proportional to each other, i.e., 
w^^^/p^^^ = w^^^/p^^^ = w^^^ / p^^^ . Therefore, w^^X'^^q^w^j)* vanishes for c = 
2, 5, 7, and 

y^i^Tci^Ziy ^ y^i^Tci^ZT ^ ^i>^Aci^Zy ^ 2UUUai)^ = 4 (3 9) 

plp2 plpi p2p3 

Furthermore, similar terms with c = 3 and c = 8 are given by 

^f.ArJ(^S)* = 2 {{p'f - {p'f) , wiX^^^i^ = -^ ((pl)2 + (p2)2 _ 2(p3)2) . 

Using the symmetric structure constant rf"^^, we get 
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(3.10) 



Next, we consider 2AqAqJ°-. Again, due to the simplicity of the ADHM data for the 
't Hooft instantons, j" is nonzero only for a = 1,4, 6, 3, 8, and we obtain 

2A^,A',r = AAlA^p' + AAtAyp' + 4A^A°pV' 

+2AIA', {{pr - ipr) + ^« iipr + ipr - ^ipr) -(s.n) 

In addition, we have {Aff (Ea \^i?) = 2 ((p')' + (p')' + (P^)') (^o)'- So, in total, 
the kinetic term for w is evaluated as 
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+p'p'AlK - ^« + P"P"« + P'p'^l^o - P'P'AIAI 

+ {-^^i^i + \i<r + \{Air - \{Air - \{Aiy) Upr - ipr) 
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X ((pY + (P'f - 2(p^n . (3.12) 

Finally, the CS term has only the A^ element, 

1627r ,n 
^cs = T71— ^0- (3.13) 

AMkk 

The total action Laq involving the gauge field Aq is a sum of ( |3.4|) , ( p.l2|) and (|3.13|) , 
as 

Lao = tT{DoX^'f + tiDowtDoWai + Lcs (3.14) 

Procedure (3): integrate out Aq 

We have to solve the simultaneous equations for all A'^ and Aq, 

^ = ^=0 (foralla = l,...,8) (3.15) 

Although Aq is mixed with the other components A^ a unique solution is found as 

° XMkk {rl^p^p^ ' ° AM^;, (r|^ + v/3rf0 W ' 

Al^l''; y^'^^P' , Al^Al^Al^O, (3.16) 

° XMKK{ri'-VSri'yp^p^ o o o > v ; 



Aq, Aq and A^ have complicated expressions, so we omit to write them exphcitly 
here. We plug the solution back to the action Laq given by ( p.l4| ), then we obtain 



integrated action Laq in terms of the moduli parameters r 
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The total Hamiltonian (potential) V^^ is given by 
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as in the two-body case. We obtain 
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(3.19) 



To find the potential intrinsic to the three-body, we need to subtract the one- 
body and two-body Hamiltonians. For the ADHM data for the 'tHooft instantons, 
they are given by |jl| 
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where r is the distance between the two baryons. The subtraction of these give 
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(3.21) 



A=l,2,3 A^B 

Therefore, the three-body forces vanish, for the baryons sharing the same classical 
spin/isospins. This result is the same as the one given in the soliton approach p]. 

The "classical" spins and isospins are realized when the magnitude of the spin/isospins 
is large, which is only possible for heavy higher spin baryons, but not for spin 1/2 
nucleons. Therefore unfortunately this "classical" treatment does not work for the 
realistic nucleons. Next, we keep the quantum spin/isospin degrees of freedom (the 
phase in w) explicitly in the computation and provide a framework for nuclear forces 
with standard quantum spin/isospins. 
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3.2 Generic three-body interactions: a set-up. 

Procedure (1): solving the ADHM constraint 

The computations with the ADHM data for the 'tHooft instantons are easy but 
they are not realistic system, since the spin/isospin rotation matrix U is fixed by 
hand. We have to allow arbitrary U for each baryon, in general. This means, instead 
of the previous (|3.1|) , we allow'^ 

y^i = U^^P^ (A = 1,2, 3). (3.22) 



In order to satisfy the ADHM constraint ( |2.2| ) with this generic w, the off-diagonal 
components of the matrices X*^ should be turned on, instead of ( p.2|) . 



X''=j:^r!:'+ E ^r^. (3.23) 



2 " ^^ 2 

a=3,8 a=l,4,6 



The diagonal r^ and rg specify the positions of the three baryons, while the off- 
diagonal ri,r4 and rg are small. 

Although generic three-instanton ADHM data is not available, we may need only 
the ADHM data for well-separated instantons. 



^3 



+ y3r8|/2, |-r3 + v^r8|/2, |r8|>p, (3.24) 



since the classical size of the instanton (baryon) is quite small as p ~ 1/vA for large 
A. The ADHM data for the well-separated instantons is described in ||T^. In our 
notation, it is^ 

|"12| 

p,P2pldfiaM r((f;3)tf/2 _ (f;2)tf;3) ^t^^^^ ((f/l)tf/3 _ (f/3)tt/l) 



4M12PM13PM23P 

- {{u^yu' - (uyu^) 4i42 {{uyu^ - (u^yu^)] + o{i/d') , (3.25) 



^Note that in single instanton case, we have gauge freedom from Aq to chose this dow = gauge. 
In three instanton case, we have Aq, Aq, and Aq gauge freedom to choose this dow^^^ — dow^^^ — 
dow^=^ = 0. 

■^Our ri/2, r4/2 and re/2 correspond to 612, &13, and &23 of [Q, as they are the off-diagonal 
elements of the matrix X^'^ . Our piU^ corresponds to qi of [Im. Our formulas (3.25), (3.26), and 



( 3.261) can be obtained explicitly from Eq. (5.13) of 1 16 , by substituting recursively the expression 



of 6,; 
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r^(^M 



'Jd^ 



p,p, {{u'^yu' - (f/i)tf/3) 



Plp3pid^3(^M 



(([/2)tf/3 _ (f/3)tf/2) 4^21 {{uyu' - (uyu') 



Md,2\'\dis\^\d23\^ 

■ {{uyu' - (uyu^) 4it/23 {{u^u^ - (uyu^)] + o{i/d') , (3.26) 



rfau = ^P2Ps iiU^U' - (UyU^) 



M2-I2 



^23 1 



p2p3pid23(^M 



41/7 |2|/7 121^ |2 
4|ttl2| |"13| |"23| 



((f/l)tf/3 _ (f/3)tf/l) 4^^^ ((f/2)t^l _ (f/l)t;72) 

- ((f/i)^f/2 - (uyu^) d\^dr3 {{wyu' - (uyu^)] + o{i/d') . (3.27) 

Here we have defined 



dij = d,, au 



(3.28) 



where (i^ is the distance vector between the i-th and the j'-th instantons. From 
( p.23| ), the location of the first, second, and third instanton is 



M 



M 



/2 + r*V2v^, -r3*V2 + r*V2v^, -rl^/V?>., 



respectively. Therefore we have 



di2 — —d2i — T3 , 

di3 = -d3i = (rs - y3r8)/2 , 

d23 = -d32 = (-rs - V3r8)/2 . 



(3.29) 



(3.30) 
(3.31) 
(3.32) 



Procedure (2): substitute the ADHM data to the action 

As all U^ matrices are different, we need to consider w^^X 



ai^^ACV'^ai 



C\* 



for all c 



1, • • • ,8. However, because U^^ G SU(2), w^i^^Ac{''^'ii)* ^^ c = 2, 5, 7 vanish as they 
are proportional to U^^iU^^Y — U^i{U^y with A,B = 1, 2, 3. The other components 
are calculated as follows: 



wi^\ciwZy 



,A \8 



C-- 



'^ai^Aci'^at 



2 {{p'r - ip'f) 

p'pHuy{uiy+uL{uyy)- 

pViuUulV + uUum- 

^ iip'r + ipr - 2ipr) ■ 



¥p'ur 



^P'p'ur 



(3.33) 
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Here we have used the fact that the SU{2) matrices U^iiU^^y for A^ B, {A, B) = 
1, 2, 3 can be written as ^012x2 + i J2i=i '^i''"* with J2i=oi'^i)'^ — l? i^ terms of Pauh 



Ui^iUf;)^ ^vi^''{l,.^^, + ^J2v!^'%. (3.34) 



The definition of Uq follows that of the two-baryon case, ( |2.10D . 



In short, compared with the previous ADHM data for the 'tHooft instantons, we 
have new parameters r^ with a = 1, 4, 6, and u^^. We first describe the integration 
of Aq, followed by the explanation of the potential due to the mass term tr[(X^)^]. 

Let us write the terms including Aq in the matrix model action explicitly. They 
are the kinetic terms of X and w, and the CS term. Note that due to the fact 
that 'wfiA^^(ti;^j)* = again for c = 2,5,7, the calculation for the w kinetic term 
Dw{Dw)* is very similar to that of the ADHM data for the 'tHooft instantons. On 
the other hand, the kinetic term for X contains terms in ( p.4[ ) as well as terms 
including r^ with a = 1, 4, 6. 



tr(Do^*')' = tr -i 



C=1.3,4,6,8 



J.M 



(3.35) 



Due to the fact that terms including A'^ with a = 2, 5, 7 decouple from the terms 
including A^ with b = 1,3, 4, 6, 8, and the fact that all Aq with a = 2, 5, 7 appear in 
the Lagrangian as quadratic terms, the equations of motion for Aq with a = 2, 5, 7 
are simply solved hj Aq = {a = 2,5,7). With this observation, the kinetic term 
for X^^ is simplified as 

tT{DoX''f 

= I [Kfrl + iKfrl + AiAlYrl + {At^frl + (A^)^^ - 2AjA^rir4 

-2V3AlAlnr, + GAlAlr^r, - 2V3A^QAlrsr, + (AlfinY + iA^)\r,y 
+3{Alfinf - 2AlA'Qrire - 2V3AtAlnre - eAjA^rgre + 2V3A^QAlrsre 
-2A^QAlnre + AVsAlAlnre 

+iAir{r,r + {Atnr.r + 3{Air{r,r + {Airmnr + {nr + {r,r) 

+4V3A'QAlnrs + 2\/3(A^)V3r8 - 2V3{Al)^rsrs - 2V3AjA^r4r8 
-GA^QAlnrs - 2V3AlA^Qr,rs - eA^A^rgrg + SiA^Q^rl + 3{Alfrl 
-2AJ^(4Ajrir3 + SA^Q^n + A^qT^u - ^3^0(^4)' 
-SA^Q^re + Alr^r^ + VsAHr^f + V^A^qT^ts - V^A^rgrg)) . (3.36) 
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The kinetic term for w is similar to the previous case with the 'tHooft instanton 
ADHM data. 

\ ??=1,3,4,6,8 / 



,2^3„,(23) .6 aO 



+4pv^«r^^^^^ + ^p'p'u^o'''K< + 2AIA', {{p^r - ipr) 



+^^x iip'r + ipr - 2ipr) 

„1 „3„,(13) 
P P % .4.8 , ^ln3,,(13) .1 .6 , n^n^.M^^A^A^ 



O 1 2 (12) 

2P P Up ,1 ,8 , 1 2 (12) .4 .6 
^= /Iq^O + P P ^0 ^0^0 



P P "o 46 



+ 






4 



+^ (2K)^ + 2(A^)^ - 2iAir - (Atr - (Air) m' + ipr - npr) (3.37) 



With the CS term given by ( p.l3| ) the total Lagrangian is again written in the form 
( |3.14| ). Again, we have to solve the simultaneous equations ( |3.15| ). 

Next, we shall describe the potential due to the mass term of the matrix model. 



547r 3 L^ ^ J 



XN.Ml 



KK 



3% 



\{rl + rt/Vsr + \{-rl + rt/V^f + ^(r^^)^ 



+ ^ {{rtr + {rir + {rtf + {rlf + {rtf + {r^f) 



. (3.38) 



The first three terms correspond to the square of the diagonal elements X| for 
2 = 1,2, 3, so these correspond to the three copies of the one-baryon potential. The 
terms in the second line are the two-body and the three-body terms. To evaluate 
these, we need explicit expressions for the off-diagonal ri, r2, r^, rs, rg and r^. 

In principle, it is a straightforward calculation to determine the three-body force 
from this. However the actual calculation turns out to be extremely messy, and it 
is hard to get a physical interpretation from that. To extract the physical essence, 
next we will choose a particular alignment of the baryons to simplify the expression, 
so that final answer is easier to analyze. 

3.3 Three-body Hamiltonian for baryons aligned on a line. 

We will find that the Hamiltonian is simplified significantly when all the baryons are 
aligned on a line. We consider the following case 



,M 



0, 



,M 



r'' ^ . 



(3.39) 



14 



This means that the first, the second and the third baryon are placed at x = r^ /2, 



X 



M 



-rg /2, and x = 0, respectively. In this case, since Vg = 0, the expression for 
the potential is significantly simplified. Still, since we are treating the spin/isospin 
moduli quantum mechanically, this gives lots of information on the three-body forces. 
In addition, we notice that all the size moduli pi can be taken to be a classical value. 
Pi = P2 = Ps = P, since we are dealing with only the leading order in the large N^ 
limit. The large Nc limit is the same as classical limit since the action (|2.1| ) has 
overall Nc- This simplifies the computation too. 

The resultant Lagrangian concerning the gauge field Aq is 



Mo 



XMkkNc 
547r 



{L, + L2] 



(3.40) 



where 



162A^7r 



(^l^2.2 



r- (A^)V2 (A6)V2 



AMi 



+ 



+ 



KK 



+ {Qi<r + (Air + {Air + {A^r + {Air + {Air) p' 






p'uT' + ( AlK + KAl 






2 (13) 
P n 



46/I8 



V^ 



2 (23) 



ThAAl [AlurhAlurhAlur^) p^ 



4L2 = 2{Al 



3\2/M\2 



1 



(^0 



A\2/^M\2 



1 



+ 2(-4o 



6\2f M\2 _ A\ Ai M M _ q 43 Abr^Vl 



4 



-v/3A^A«rf rf + i(A^)^(rf )^+ ^(A^)^(rf )^+ ^(A^)^(rf )^+ VSA^A^lrf )^ 
3 

■2' 



-:^(A^)^(rf)^ + 3AgA^rrrf 



4I 46 M M _ ./q 4448 71 



M„M 



r, 



6 



/14/.6 M M 



-2V3AlAlrrrr + l{Al 
3 



-VSA^A 



8f^M\2 



o^^ov'^e 



:K 



^^^^*^)^ + ^(^^)^(rf)^ + i(^^)^(r6*0^ 



>/\2 



(3.41) 



For getting this expression of Laq = Li + L2, we have used the equations 



'3 '1 







,A/ 



+ V^r*>f = 0, (r^ 



3 



v^r^^Orf 



0, 



(3.42) 



to eliminate cross terms between rs § and y in the Lagrangian. These can be shown 
explicitly using the solution of the ADHM constraint ( |2.2|) in the expansion of the 
small p^ jr^ . The expansion was studied in detail in [|16|. Using the expression given 
in Eq. (5.13) of [|T6|, it is easy to show the equations above. This elimination of the 
cross terms is important for a simplification of the computations. In fact, we can 
show later that L^ is not necessary, when integrating out A^. 

Next we evaluate the contribution from the mass term (X^)^. The alignment 
( ^^391) simphfies the ADHM data (jpl), (13^261) and ( ^^2^ ) quite a lot. In fact, we 
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have 



'^f ^M = r^PiP2 rT2i - t^PiP2pI r{T^2Ti:, - ^^^2) + Oil/\rf) , (3.43) 



<Vm = T^PiPs rTgi - T^PiPspi r(T32Ti2 - T12T32) + 0(l/|rn , (3.44) 



r^^M = -n^P2P3 rT32 + rilP2P3Pl r{TsiT2i - T^iT^i) + 0{l/\r\'>) , (3.45) 



where r = t^Vm, and T^j = (f/*)tf/J' - (t/^)tf/^ = -Tji. Note that the first term 
in each of the right hand side equals the off-diagonal entry of the two-body case, 
Y in ( |2.7|) . The second terms are corrections due to the three-body effect. So, the 
three-body contribution in the mass term tr(X'^)^ (|3.38| ) should appear at the leading 
order as a liner term in these second terms, multiplied by the first terms. An explicit 
computation leads to 

^ma.s _ AAT.MIk p' 



3-body - 223% |r|6 

X (tr[rT2i]tr[r(r23ri3 - T13T23)] - 2tr[rT3i]tr[r(T32ri2 - T12T32)] 

-2tr[rT32]tr[r(T3iT2i - T21T31)]) . (3.46) 

We have already subtracted the one-body and the two-body potentials here, and 
took pi = P2 = P3 = P which is the classical value (the leading value in the large Nc 
expansion) . 

Procedure (3): integrating out Aq 



Once we solve the simultaneous equations of motion for ( p.l5|) for Aq, and plug 
the solutions into the Lagrangian Laq, we should obtain 

Lao = -V, V^ Y. ^l-body + ^ E ^2^oi + ^3-body (3.47) 

A=l,2,3 Aj^B 

where the first term is the one-body rest energy, and the second term is the two-body 
interaction potential. As obtained in [H], their expressions are 

y(A) 277rAr, 1 (^,^) 27vrAr, {ui''^^f 

l-body 4AMKKP2' 2-body AMkk |r(^B) |2 + 2p2 _ 2(4^^^) V 

Here the inter-nucleon distance is, according to our alignment ( p.39| ), 

|r(i2)|=^, |r(i3)| = |^{23)|^^/2. (3.49) 

The third term V3_body is what we like to compute. 

We are interested in the regime of short distances r ^ 1/Mkk- However, as the 
classical size of the baryon p is quite small, p ~ 1/(vAMkk), the region of our interest 



16 



is rather a "long-distance" expansion p ^ r in effect. Therefore we need to expand 
the resultant Hamiltonian for small p/r. As we look at the Lagrangian L^^ = L1+L2, 
we notice that L2 is of order p^ jr^ ^ as we know that the ADHM constraint is solved 
in this expansion as |/ = O^p^ /r). On the other hand, As is obvious from (|3.47|) , 
the two-body interaction is 0(l/r^) so the three-body interaction should start from 
p^ jr'^. (This is suggested also from the sohton approach, see p.) Therefore, L2 is 
not necessary as it is at higher order in this expansion.^ 
The Lagrangian L\ can be conveniently written as 



U 



P^MA^B^A 



(3.50) 



where 



TtT 



AM, 



A^ ^ ( A^ A^ 4^ 44 A^ A^\ W 

2 
M = P + Q, P= ^diag(0,4,0,l,l,0) 



1627r 



KK 



(1,0,0,0,0,0), (3.51) 

(3.52) 



/ 6 



Q^P' 



2m, 



(12) 



2u, 
2m, 



(13) 



(23) 



V 



2m, 



(12) 








1 

MfV2 4^^V2 
4"V2 -nrV2 



u 



(12) 



/Vs 







2ur' 

MfV2 
4"V2 

4"V2 

,(13) 



-nrv2v^ -«rv2v^ 



2Mf) 

4"V2 

-4"V2 

4"V2 

1 

(23) 



U 



-u, 



-u. 



fVv/3 



\ 



/2x/3 



(13) 

fV2v^ 
1 



(3.53) 



Since M is a symmetric matrix i.e. M^ = M, the equations of motion for Aq is 
solved by 



A = — M-^B, 
2 

which is substituted back to Li to give the Hamiltonian (which is —La^] 

547r 4 



(3.54) 



(3.55) 



As B has only one non-zero entry, this is nothing but 



V 



S^irNr 



2AMi 



KK 



[M-1 



J (1,1) 



(3.56) 



^As a check, we can perforin a computation with keeping L2 exphcitly to confirm this. The 
computation is lengthy and is not presented in this manuscript, but we have confirmed it. Note 
that for the generic case with nonzero rg, this simphfication is not expected, because in general there 
are terms of the form r^y which contributes additionally to Li, so one needs explicit expression for 
y by solving the ADHM constraints. 
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which can be evaluated using the first cofactor of the matrix M. By expanding in 
power series of p^ /r"^ up to 0{p*/r^), we obtain, 



where spin/isospin phase /si is defined as 

-244^'^)n?'^)4^'^) . (3.58) 



Subtracting the 1-body and 2-body potentials (|3.48D from this expression as in (|3.47|) 



2 /^2 



we obtain the potential intrinsic to the three-body nature by the expansion of p /r 
which we call V^.^ody ^^ 



K 



^„ _2167riV^ 



3-body AMKK|r|4 



rf'^Vrf'^V + rf'^y rf'^y + 4(4^'^V(u?'^V 



-64^'^)t.(^'^)4^'^) 



+ 0{pyr^) . (3.59) 



With the potential coming from the X^ mass term (p.46 ), the total three-body 
potential is 

V's-body = V'g.body + ^3-body • (3.60) 

With this at hand, we can evaluate this potential with any three-baryon state with 
any spin/isospin. Next, we shall choose two wave functions, one is appropriate for 
the neutron stars, and the other is for a Helium-3 nucleus and a triton (a nucleus of 
tritium), to find the three-body nuclear potential is positive. 

We have two remarks on ( [3.591 ). First, the three-body Hamiltonian ( ^.59] ) is of 
order (9(l/(A^r^)) because of p ~ 1/vA, and so it is suppressed by 1/A^. This is 
consistent with the generic observation given in the soliton picture |^ stating that the 
generic fc-body potential is of order l/(A'^~^r^'^~^) in the unit of Mkk = 1- Second, 
in the expression above if we put all the matrices f/*^*-* equal to each other so that 
the ADHM data is that of the 'tHooft instantons, we have Uq = 1 and Aij = 0, 
resulting in the vanishing three-body potential. This is consistent with the result of 
the previous section. 

Procedure (4): evaluate the Hamiltonian with baryon states 

Now, we are ready to compute the spin/isospin dependence of the three-body 
nuclear force at short distances. Although we can evaluate it for any choice of 
spin/isospin for each baryon, in this paper we choose the following two states as 
exphcit examples: 
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(4-a) three neutrons with spins averaged. 

(4-b) proton-proton- neutron (and proton-neutron- neutron) . 

The reason for these choices is that the first example is relevant for dense states of 
many neutrons such as core of neutron stars and supernovae, where the three-body 
nuclear forces are quite important. The second is obviously for the spectrum of 
Helium-3 nucleus where three-body forces are expected to contribute, and also for a 
triton. 

(4-a): three neutrons with spins averaged 

For protons and neutrons, the single-baryon wave function is given by (|2.5|) . For 
the neutron stars and the supernovae, we need neutron states with spins averaged. 
Thus, for any given operator of the quantum mechanics, the appropriate expectation 
value for these is obtained by 



{V) = - {n\\0\n\) + {ni\0\n\) 



(3.61) 



For the case of O being the three-body Hamiltonian, we need to take the above 
expectation value for each of three baryons. As nucleons are fermions, any wave 
functions should be anti-symmetric under the exchange of the nucleons. Here, as 
three neutrons move around in realistic situations, we do not anti-symmetrize the 
wave functions^ (in this paper we have not evaluated nuclear potentials coming from 
motion of the baryons) . 

Here for a demonstration, let us consider a single baryon case ( p. 611 ). Using the 
coordinate expression of the wave functions, this ( |3.61| ) means 



{V) 



da, - 



0\{d\n^)\^ + 0\{d\nl)\' 



(3.62) 



Here dfl^ is the integration over the S^ spanned by the unit vector a. Using the wave 
functions (|2.5|) , we can see 



(a|nt)r + |(a|^4)P 



TT^ 



[ai? + (a2? + (as? + (a^? 



1 



TT' 



(3.63) 



So, we obtain a simple expression 

{V) 



27r2 



dQ^O. 



(3.64) 



^In fact, once we take three neutrons for the isospin sector, it is impossible to anti-symmetrize 
the wave function with the spin sector, without a help of angular momenta. 
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Using this simple formula, the three-body potential with the spin-averaged wave 
function is 



(•^S-body/ 



body / nnn(spin— averaged) 



{uS'yiuS''r+iuS''nu^^''r 



AMKK|r|4 (27r2)3 J -^"^ ^ ^ 

+4(4^'^V(4^'^V-64^-^)4''^)4''^)l. (3.65) 



This integral over three S'^'s can be easily performed. For example, for (-Ug ' )^, using 



(1,2)n2 



the definition below ( |2.5|) and ( p.lO|) , we get 



M^''^') = -tr [U^'^^U^^^] = a« ■ a(^') , (3.66) 

where a'-*-' is unit 4-component vector, pointing one phase point on S^ for spin/isospin 
d.o.f. U^^\ Therefore, we obtain 

fdni^\u^Q'^^f = fdn'i^ cos^ 9 = /"cos^ 9 sin^ 9 sin 9d9d9d9 = — , (3.67) 
where 9 is the angle between a^^^ and a^'^^ . Using this and also the following integral 

fdQ'i^dQf^dQ'iK^^'\f''\^^''^ = y , (3.68) 

we obtain 

\ J— body /Qmj^gpjjj_^ygj.3^g(;(j^ v / 

Therefore the three-body potential from the Aq term vanishes for the spin-averaged 
neutron wave function. 

In the same manner, for V^^b^ody ^^^ expectation value is given as 



/t Amass \ 
Vl^a-body/ 



nnn(spin— averaged) 2^3^77 Irl 



X j^ Jdn'i^dnfdnf'> [tilrT^iMriT^sn, - T,^T^,)] 

-2tr[rT3i]tr[r(T32Ti2 - T12T32)] - 2tr[rT32]tr[r(T3iT2i - T21T31)] 



22337r(27r2)3 Irl 



jdn^^UnfUnf (tr[rT2i]tr[r(r23Ti3 - T13T23)]) . (3.70) 



Here in the last equality we have used the invariance under the exchange of the 
integration variables, dVt^, O dVlr^ <H- dVt\^ . The integration can be performed by 
using the polar coordinates of the S*^, and the result is 

-i— jdn^^Unfdnf (tr[rT2i]tr[r(T23Ti3 - T13T23)]) = -8|rf , (3.71) 
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where r = {r^,r'^, r^) is the three-dimensional vector which specifies the inter-baryon 
distance in our space. At the leading order in 1/A^ expansion, we may use the 
classical value for r^ which is zero, so in eff'ect the three-dimensional distance is 
equal to the four-dimensional one, \f\ = \r\. We can substituting the classical value 
p = 2~^/^3^/^A/7rA~^/^Mj^j^ at the leading order in the 1/Nc expansion. So, we obtain 
the three-body potential due to the matrix model mass term as 

9-1/2315/2 2/Y 
\ 3-body/nnn(spin-averaged) A^M^ Irl^ ' V"-"- ' ^7 

Therefore, in total, we obtain 

\ 3-body/nnn(spin_averaged) \ 3-body/nnn(spin-avcraged) \ 3-body/nnn(spin-averaged) 

2-l/23l5/2^2jy 



A2M3, 



r 



4 



(3.73) 



'KKI 

This is the three-body nuclear potential for three neutrons placed on a line with 
equal spacings |r|/2, with spins averaged. The three-body potential is suppressed 
compared to the two-body potential by ~ l/XlrM^x)^ ^ 1 for large A, which is 
generic hierarchy between A^ + 1-body potential to A^-body one as shown in [^ . Mkk 
roughly indicates the QCD scale, and our computation is valid at short-distance, 
1/(\/AMkk) < \r\ < 1/Mkk.' 

(4-b): proton-proton-neutron 

Let us evaluate the three-body potential with the case of proton-proton-neutron. 
We are interested in the three- nucleon state with a total spin 1/2 and a total isospin 
1/2. For any choice of the third component of the spin/isospins, we can find a unique 
wave function with a complete anti-symmetrization. 

The proton-proton-nuetron means the third component of the total isospin is 
-1-1/2. For example, when the third component of the total spin is -1-1/2, 

-7^ \P t)lb l)2\n t)3 - \P i)l\p t)2\n t)3 - \P t)l|^ t)2b i)3 

+ \P i)i\n t)2b t)3 - 1^ t)ib i)2b t)3 + \n t)ib t)2b Ds] ■ (3.74) 
The calculation with this wave function is straightforward, and we find the integrals 

jdQfUnfUnf'^lijia,, d2, a,)\\'i''^u^i''^u'i'''^ = ^ , (3.75) 

y"rfn(^)rfnf)rffif |^(ai,a2,a3)p(tr[rT2i]tr[r(T23Ti3 -T13T23)]) = -^kT- 



^This Mkk is about 1 GeV if it is fit with p meson mass Q, while it is about 0.5 GeV when 
it is fit with baryon mass differences [^. We are working in the large A expansion. The 'tHooft 
coupling constant of QCD, A, is C(fO — 20) when it is fit with pion decay constant jg). 
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Using these formula, we obtain again 

(^3lU>pp„ = 0' (3-76) 

while for the other potential we have a different factor 

95/239/25 2 /Y 

\^3-body;pp, - A2M3j^|r|4 ' ^ > 

Therefore, in total, the three-body potential is 

26/239/25^2^^ 

(V3-body)ppn = C^S-body/ppn + \^3™body)ppn = _)^2j\^3 1^14 " (^•'^^) 

The three-body potential is positive, that means, we have a repulsive three-body 
force at short distances. 

This computation is for (+1/2, +1/2) of the third components of the spin and the 
isospin. Computations with three other wave functions, (+1/2, —1/2), (—1/2, +1/2), 
and (—1/2, —1/2), can be done in the same manner, and the result for the three-body 
potential turns out to be the same as ( p.78| ) for all of these. These are due to the fact 



that the action (|2.1|) has rotational invariance 5*0 (3) and isospin SU{2) invariance. 
This includes the case for proton-neutron-neutron, which is the case for a triton (a 
tritium nucleus). 

4. Summary and discussions. 

With the simple U{k) matrix model for fc-nucleon systems which we proposed in 
[0] together with P. Yi, in this paper we have computed short-distance three-body 
nuclear forces. Our matrix model is not a phenomenological model, but derived in 
string theory using the gauge/string duality (the AdS/CFT correspondence). More 
precisely, our matrix model is a low-energy effective field theory on baryon vertex D4- 
branes in the D4-D8 holographic model |^, |^ of large N^ QCD. In this framework, 
we can compute nuclear potentials for arbitrary number k of the nucleons. 

Our computations are straightforward. For three nucleons, we took k = 3, 
i.e. U{3) matrix model. The matrix model Hamiltonian evaluated with a quantum 
three-baryon state, a tensor product of single-baryon states, gives the three-body 
nuclear potential. We subtracted one-body and two-body contributions, thus the 
remaining is the force intrinsic to the three-body. The computations are valid only 
at short range, l/y/XM^K *^ kl ^ 1/^kk where A is the 'tHooft coupling constant 
of the QCD which is 0(10) for fitting pion decay constant [|], and Mkk is 0{1 — 0.5) 
GeV [§, [1^ when it is fit with meson/baryon masses (or mass differences). As 
explicit examples, we took a) three neutrons with spins averaged, and b) proton- 
proton-neutron, both aligned on a line with equal spacings. The resultant three-body 
nuclear potentials are ( p.73|) and ( p. 78]) , both of which are positive. 
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Let us discuss possible importance of our result. We have computed ( |3.73| ) for 
three-body neutrons. But as seen from the form of wave functions ( |2.5| ) and isospin 
SU{2) invariance of the action ( |2.1| ), the results hold also for three-body forces for 
three-protons. Therefore the result ( |3.73| ) hold as far as all three nucleons have same 
flavor. In the same manner, the three-body potential (|3.78|) for the proton-proton- 



neutron is equal to the three-body potential for proton-neutron-neutron, which is 
responsible for a triton. These results imply that there are additional repulsive 
forces in addition to two-body forces for these states at short distances. 

The three-body potentials which we obtained in ( 3.73|) and ( |3.78 ) are suppressed 



by l/A(rMx_ft')^ compared with two-body potential, and at the length scale where 
our computation is valid, i.e. 1/(vAMkk) ^ |^| ^ 1/^kk, this suppression fac- 
tor 1/A (rMi^x)^ is small. This makes our three-body potential computation valid; 
we have two-body dominant repulsive potential and furthermore small but nonzero 
repulsive potential from three-body forces. 

These three-body forces are stronger as distances get shorter. As a result, at very 
short distances where neutrons are highly dense, three-body forces give additional 
repulsive forces. This statement supports recent observation that the nuclear two- 
body repulsion is not enough to explain supernovae explosions, nor the equations 
of states for the core of neutron stars. In high density nuclear matters such as the 
neutron stars, our result suggests that the repulsive core of neutrons in neutron stars 
and supernovae has an extra positive contribution besides the repulsive potential 
from the two-body nuclear potential. The necessity of the repulsive three-body forces 
for neutrons has been indicated by analysis of mass bounds of neutron stars and 
supernova explosion simulations. 

We also found that the three-body forces for proton-proton-neutron at short 
distances is repulsive. In IIelium-3 nuclear spectrum, it is expected that a short- 
range repulsive three-body forces is necessary, and our result sounds to be consistent 
with this. Furthermore, we found a repulsive three-body forces also for a proton- 
neutron-neutron, which should be related to a triton. There are other related issues 
in few-body nuclear spectra*^. 

Our example is limited to three nucleons on a line, so this is not conclusive for 
the questions concerning the interesting situations listed above. However, our results 
are suggestive. In principle, it is very straightforward to compute the A;-body forces 
at arbitrary arrangements of nucleons using our matrix model, therefore this matrix 
model is effective for studying short-range many-body nuclear forces. 
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